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Sensitivity of KASCADE-Grande data to hadronic interaction models 
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KASCADE-Grande is a large detector array dedicated for studies of high-energy cosmic rays in the 
primary energy range from 100 TeV to 1 EeV. The multi-detector concept of the experimental set-up 
offers the possibility to measure simultaneously various observables related to the electromagnetic, 
muonic, and hadronic air shower components. The experimental data are compared to predictions of 
CORSIKA simulations using high-energy hadronic interaction models (e.g. QGSJET or EPOS), as 
well as low-energy interaction models (e.g. FLUKA or GHEISHA). This contribution will summarize 
the results of such investigations. In particular, the validity of the new EPOS version 1.99 for EAS 
with energy around 100 PeV will be discussed. 



I. INTRODUCTION 

Studies of high-energy cosmic radiation by means 
of extensive air shower (EAS) techniques require a 
proper understanding of high-energy interactions in 
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the Earth's atmosphere: Inferring the properties of 
the primary particles from EAS measurements one 
relies on the simulations of the air shower develop- 
ment, whose backbone is the hadronic cascade in the 
atmosphere. Significant progress has been made dur- 
ing recent years to interpret air shower data and main 
physical properties of the primary cosmic ray parti- 
cles have been measured. In the energy range around 
10 15 eV, energy spectra for elemental groups and mass 
compositions of primary cosmic rays have been inves- 
tigated. Interpretations of these extensive air shower 
measurements are generally related to air shower mod- 
els to obtain physical properties of the shower includ- 
ing primary particles. Therefore, one of the goals of 
KASCADE-Grande is to investigate high-energy in- 
teractions in the atmosphere and to improve contem- 
porary models to describe such processes. 

The tests of hadronic interaction models require 
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detailed measurements of several shower observables. 
The KASCADE-Grande experiment with its multi- 
detector concept of the experimental set-up, measur- 
ing simultaneously the electromagnetic, muonic, and 
the hadronic shower components, is particularly de- 
signed for such investigations. 

The KASCADE array measures an extensive air 
shower in the energy range of 10 14 to 8xl0 16 eV and 
consists of 252 scintillator detector stations with un- 
shielded and shielded detectors located on a grid of 
200 x 200 m 2 for the measurement of the electromag- 
netic and muonic shower components independently 
In its center, an iron sampling calorimeter of 16x20 
m 2 area detects hadronic particles. 

The KASCADE-Grande [1] array covering an area 
of 700 x 700 m 2 is optimized to measure extensive air 
showers up to primary energies of 1 EeV. It comprises 
37 scintillation detector stations located on a hexago- 
nal grid with an average spacing of 137 m for the mea- 
surements of the charged shower component. Each of 
the detector stations is equipped with plastic scintil- 
lators covering a total area of 10 m 2 . 



II. MONTE-CARLO SIMULATIONS 



The principal idea of the tests of hadronic inter- 
action models is to simulate air showers initiated by 
protons and irons nuclei as the two extremes of pos- 
sible primary particles and to compare those simu- 
lated showers with the measurements. For the air 
shower simulations the program CORSIKA @ has 
been used, applying different embedded hadronic in- 
teraction models. To determine the signals in the indi- 
vidual detectors, all secondary particles at the ground 
level are passed through a detector simulation pro- 
gram using GEANT package. The predicted observ- 
ables at ground level, such as e.g. the number of elec- 
trons, muons and hadrons are then compared to the 
measurements. 

FLUKA || and GHEISHA [4j (E < 200 GeV and 80 
GeV, respectively) models have been used for hadronic 
interactions at low energies. High-energy interactions 
were treated with different models QGSJET-II-2 [|] 
and EPOS 1.99 [6j. Showers initiated by primary pro- 
tons and iron nuclei have been simulated. The simula- 
tions covered the energy range of 10 14 - 10 18 eV with 
zenith angles in the interval 0° - 42°. The spectral 
index in the simulations was -2 and for the analysis 
it is converted to a slope of -3. The simulated events 
are analyzed by the same method as the experimental 
data, in order to avoid biases by pattern recognition 
and reconstruction algorithms. 



III. INVESTIGATIONS WITH KASCADE 

Using KASCADE measurements, the hadronic in- 
teraction models QGSJET @] (versions 98 and 01) , 
SIBYLL EM (versions 1.6 and 2.1), DPMJET [lj|, 
VENUS [ljjand NEXUS [3 have been investigated. 
First tests [13[ supported QGSJET 98 as being most 
compatible with the KASCADE data. Similar conclu- 
sions have been later drawn for QGSJET 01 14]. For 
the next model version, QGSJET-II-2 [5j, some prob- 
lems with the electron-hadron correlations have been 
observed [la ]. Predictions of SIBYLL 1.6 were not 
compatible with air shower data, in particular there 
were strong inconsistencies for hadron-muon corre- 
lations. These observations improved the develop- 
ment of SIBYLL 2.1 19! . The predictions of this 
model are very successful and fully compatible with 
KASCADE air shower data [T3|. Investigations of 
the DPMJET version 2.5 yield significant problems 
particularly for hadron-muon correlations, while the 
newer version DPMJET 2.55 is found to be com- 
patible with air shower data 14j . The predictions 



of the VENUS model revealed some inconsistencies 
in hadron-electron correlations 13]. The predictions 
of NEXUS 2 were found to be inconsistent with the 
KASCADE data in particular for the investigations of 
hadron-electron correlations [3] . Presently, the most 
compatible predictions are obtained from the models 
QGSJET 01 and SIBYLL 2.1. 

Predictions of the interaction model EPOS 1.61 
have been recently compared to KASCADE air shower 
data [16]. This model is a recent development, histori- 
cally emerging from the VENUS and NEXUS models. 
This analysis indicates that the EPOS 1.61 delivers 
not enough hadronic energy to the observation level, 
and also the energy per hadron seems to be too small. 
Presumably, the inconsistency of the EPOS predic- 
tions with the KASCADE measurements is caused by 
too high inelastic cross sections for hadronic interac- 
tions implemented in the EPOS model. To solve these 
problems, the treatment of screening effects in nuclear 
collisions has been improved in EPOS. The new ver- 
sion EPOS 1.99 [(| has a reduced cross section and in- 
elasticity, compared to the previous EPOS 1.61 which 
leads to deeper shower development. 



IV. HADRONIC MODELS QGSJET-II-2 AND 
EPOS 1.99 



QGSJET is essentially based on the Quark-Gluon- 
String model approach to high-energy hadronic inter- 
actions, including a generalization of the latter for 
nucleus-nucleus collisions and a treatment of semi- 
hard processes, using the so-called semihard Pomeron 
approach |7|. The new hadronic interaction model 
QGSJET-II accounts for non-linear interaction effects, 
which allows one to obtain a consistent description of 
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FIG. 1: Two-dimensional shower size spectrum measured 
by KASCADE-Grande, together with proton and iron in- 
duced showers for QGSJET-II-2 and EPOS simulations. 



hadronic cross sections and parton distribution func- 
tions EPOS @ is a consistent quantum mechani- 
cal multiple scattering approach based on partons and 
strings, where cross sections and the particle produc- 
tion are calculated consistently, taking into account 
energy conservation in both cases. 

Predictions of QGSJET-II-2 and EPOS 1.99 models 
have been investigated with KASCADE-Grande data. 
Figure [T] represents the measured two-dimensional 
shower size spectrum, color coded area, which includes 
full detector response by simulations. The symbols 
correspond to the primary protons and iron nuclei, as 
predicted by the interaction models EPOS 1.99 and 
QGSJET-II-2. The errors of mean values are plot- 
ted in Fig. [TJ It is shown that the most probable 
values for EPOS arc shifted toward the higher muon 
numbers with respect to QGSJET. This behavior im- 
plies, if EPOS predictions are used to derive the mass 
of primary particles from the observed data, a domi- 
nantly light mass composition. Air showers simulated 
with EPOS 1.99 have about 10% more charged parti- 
cles and about 15% less muons than QGSJET-II-2 at 
KASCADE-Grande energies. 

The influences of the different hadronic interac- 
tion models on the reconstructed all-particle energy 
spectrum was investigated by performing the recon- 
structed charged particle shower size method, based 
on simulations with the hadronic interaction models 
QGSJET-II-2 and EPOS 1.99. The shower size per 
individual event is corrected for attenuations in the 
atmosphere by the constant intensity cut method and 
calibrated by Monte-Carlo simulations under the as- 
sumption of a dependence IgE oc lgN c h and a partic- 
ular primary composition. To determine the energy 
conversion relation between the number of charged 
particles and the primary energy, the simulations were 
used. The relation of the primary energy as a function 
of the number of charged particles is shown in Fig. [2] 
for the assumption of primary protons and iron, re- 
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FIG. 2: The primary energy as a function of the num- 
ber of charged particles for assumed pure proton and iron 
components for QGSJET-II-2 (top panel) and EPOS 1.99 
(bottom panel), respectively. The lines show the applied 
fits to the points. 



spectively. Assuming a linear dependence in logarith- 
mic scale of IgE = a + b-lgN c h, the fit is applied in the 
range of full trigger and reconstruction efficiencies. A 
detailed process for the energy reconstruction is given 
in Ref. 0. 

Figure [3] shows the all-particle energy spectra ob- 
tained after applying the energy reconstruction func- 
tions as well as the appropriate correction for the 
bin to bin fluctuations, based on the assumption of 
iron and proton for QGSJET-II-2 and EPOS mod- 
els. EPOS results lead to significantly higher flux (10- 
15%) compared to QGSJET. Since the EPOS model 
has for a fixed primary energy less charged particles, 
it assigns higher flux. As the calibration depends on 
simulations, other interaction models would change 
the interpretation of KASCADE-Grande data, so that 
more investigations by applying and comparing vari- 
ous interaction models are needed. 



V. MUON DENSITY INVESTIGATIONS 

The muon density can be directly measured by 
KASCADE, so that the composition studies as well 
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FIG. 3: Reconstructed all-particle energy spectra from 
KASCADE-Grande shower size for assuming proton and 
iron composition, based on two different hadronic interac- 
tion models QGSJET-II-2 and EPOS 1.99. 



as the tests of hadronic interaction models with muon 
densities can be performed . Figure 2] shows the 
correlation of muon density with the electron numbers 
and the distance from the shower axis, compared to 
the predictions of QGSJET-II-2 and EPOS 1.99 us- 
ing proton and iron nuclei as primary particles. The 
muon density decreases with increasing distance from 
the shower axis and increases with increasing electron 
number. Since an equal probability trigger for pro- 
tons and iron primaries as a function of distance from 
the shower axis is assumed, one should expect the lat- 
eral density distribution to be parallel to pure com- 
position primaries. It shows that the lateral distri- 
butions of simulated proton and iron shower are par- 
allel, while the measured one is not quite parallel to 
the both QGSJET-II-2 and EPOS 1.99 curves. How- 
ever, the lateral muon density distribution has a better 
slope than the other models such as EPOS 1.61. The 
QGSJET-II-2 model could fit the data with an inter- 
mediate primary abundance between proton and iron 
nuclei, whereas EPOS 1.99 would require abundance 
of light primary particles in order to fit the data. 



1.99 model indicates that light abundances of primary 
cosmic ray particles would be needed to fit the data. 
On the other hand, the QGSJET-II-2 model describes 
the data with an intermediate primary abundance 
between proton and iron nuclei. The reconstructed 
all-particle energy spectra are presented by using the 
hadronic interaction models QGSJET-II-2 and EPOS 
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FIG. 4: Muon density as a function of the total number 
of electrons (top panel) and lateral distribution of muons 
(bottom panel) compared to the predictions of QGSJET- 
II-2 and EPOS 1.99. 



VI. CONCLUSION 



Testing of hadronic interaction models QGSJET-II- 
2 and EPOS 1.99 implemented in the CORSIKA pro- 
gram have been performed with KASCADE-Grande 
air shower data in the energy range of 10 16 to 10 18 
eV. From the muon density investigations, the EPOS 



1.99. The resulting spectra show that the interpreta- 
tion of the KASCADE-Grande data with EPOS 1.99 
leads to significantly higher flux as compared to the 
QGSJET-II-2 result. Mored detailed investigations of 
EPOS 1.99 is still in work. 
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